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We calculate the resistance between two arbitrary grid points of several infinite lattice structures of 
resistors by using lattice Green's functions. The resistance for d dimensional hypercubic, rectangular, 
triangular and honeycomb lattices of resistors is discussed in detail. We give recurrence formulas for 
the resistance between arbitrary lattice points of the square lattice. For large separation between 
nodes we calculate the asymptotic form of the resistance for a square lattice and the finite limiting 
value of the resistance for a simple cubic lattice. We point out the relation between the resistance 
of the lattice and the van Hove singularity of the tight-binding Hamiltonian. Our Green's function 
method can be applied in a straightforward manner to other types of lattice structures and can be 
useful didactically for introducing many concepts used in condensed matter physics. 

I. INTRODUCTION 

It is an old question to find the resistance between two adjacent grid points of an infinite square lattice in which 
all the edges represent identical resistances R-r-Xt is well known that the result is R/2, and an elegant and elementary 
solution of the problem is, given by AitchisonEJ. The electric-circuit theory is discussed in detail in a classic text by 
van dejt Pol and Bremmeia and they derive the resistance between nearby points on the square lattice. In Doyle's and 
Snell'sfj book the connection between random iMa-lks and electric networks is presented, including many interesting 
results and useful references. Recently VenezianQ and Atkinson et. al.H also studied the problem and in these papers 
the reader can find additional references. Venezian's method for finding the resistance between two arbitrary grid 
points of an infinite square lattice is based on the principle of the superposition of current distributions. This method 
was further utilized by Atkinson et. al.El and applied to two dimensional infinite triangular and hexagonal lattices as 
well as infinite cubic and hypercubic lattices. In this paper we present an alternative approach using lattice Green's 
functions. This Green's function method may have several advantages: (i) It can be used straightforwardly for more 
complicated lattice structures such as body and face centered cubic lattices, (ii) The results derived by this method 
reflect the symmetry of the lattice structures although they may not be suitable for numerical purposes. However, 
from these results one may derive other integral representations of the resistance between two nodes, which can be 
used for evaluating the integrals either algebraically or numerically with high precision. Later in this paper we shall 
give examples for this procedure, (iii) From the equation for the Green's function one can, in principle, derive some 
so-called recurrence formulas for the resistances between arbitrary grid points of an infinite lattice. In this paper such 
recurrence formulas are derived for an infinite square latticeo (for the first time to the best of the author's knowledge), 
(iv) In condensed matter physics the application of the lattice Green's function has become a very efficient tool. The 
analytical behavior of the lattice Green's function has been extensively studied over the past three decades for several 
lattice structures. We make use of the knowledge of this analytical behavior for an infinite cubic lattice and give 
the asymptotic value of the resistance as the separation of the two nodes tends to infinity. Throughout this paper 
we shall utilize the results known in the literature about the lattice Green's function and also give some important 
references, (v) Finally, our approach for networks of resistors may serve as a didactically good example for introducing 
the Green's function method as well as many basic concepts such as the Brillouin zone (BZ) used in solid state physics. 
We therefore feel that our Green's function method is of some physical interest. 

The application of the Green's function proved to be a very effective [-method for studying the transport in inho- 
mogeneous conductors and it has been used successfully by Kirkpatricku for percolating networks of resistors. Our 
approach for obtaining the resistance of a lattice of resistors is closely related to that used by Kirkpatrick. Economou's 
boolfl gives an. excellent introduction to the Green's function. A review of the lattice Green's function is given, by. 
Katsura et. al.D The lattice Green's function is also utilized in the theory of the Kosterlitz-Thouless-Berezinski£3~tL3 
phase transition of the screening of topological defects (vortices) . A review of the latter problem is given by Chaikin 
et. al. in their booklia The phase transition in classical two-dimensional lattice Coulomb gases has been studied by 
Lee et. al.c3 also by using the lattice Green's function in their Monte Carlo simulations. The above examples of the 
applications of the lattice Green's function are just a selection of many problems known in the literature. A review 
of the Green's function of the so-called tight-binding Hamiltonian (TBH) used for describing the electronic band 
structures of crystal lattices is presented in Economou's booku. The lattice Green's function defined in this paper is 
related to the Green's function of the TBH. Below we shall point out that the resistance in a given lattice of resistors 
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is related to the Green's function of the TBH at the energy.-at.-»jhich the density of states is singular. This singularity 
is one of the van Hove singularities of the density of statesliScZI. 



II. HYPERCUBIC LATTICE 



Consider a d dimensional lattice which consists of all lattice points specified by position vectors r given in the form 

r = Ziai+; 2 a 2 H h h&d, (1) 

where ai, a 2 , • • • , are independent primitive translation vectors, and l\, I2, • ■ ■ , Id range through all integer values 
(i.e. positive and negative integers, as well as zero). In the case of a d dimensional hypercubc all the primitive 
translation vectors have the same magnitude, say a, i.e., |ai| = |a 2 | = • • • = |a f /| = a. Here a is the lattice constant 
of the d dimensional hypercube. 

In the network of resistors we assume that the resistance of all the edges of the hypercube is the same, say R. We 
wish to find the resistance between the origin and a given lattice point ro of the infinite hypercube. We denote the 
current that can enter at lattice point r by I(r). However, for measuring the resistance between two sites the current 
has to be zero at all other sites. Similarly, the potential at site r will be denoted by V(r). Then, at site r, according 
to Ohm's and Kirchhoff's laws, we may write 

I(r)R = ]T [V(r) - V(r + n)] , (2) 
11 

where the n are the vectors from site r to its nearest neighbors (n = ±aj, i — 1, • • • , d). The right hand side of Eq. (|^) 
may be expressed by the so-called lattice Laplacian defined on the hypercubic lattice: 

A (r) /(r)=£[/(r + n)-/(r)]. (3) 

n 

The above defined lattice Laplacian corresponds to the finite-difference representation of the Laplace operator. The 
lattice Laplacian l/a 2 A( r )/(r) yields the correct form of the Laplacian in the continuum limit i.±~ a — > 0. The lattice 
Laplacian is widely used to solve partial differential equations with the finite-difference methodEa. 

To measure the resistance between the origin and an arbitrary lattice point ro we assume that a current I enters 
at the origin and exits at lattice point ro. Therefore, the current is zero at all the lattice points except for r = and 
ro, where it is I and —I, respectively. Thus, Eq. (0), with the lattice Laplacian, can be rewritten as 



where current at lattice point r is 



The resistance between the origin and ro is 



A (r) V(v) = -I(t)R, (4) 

I(r) =J(a r ,o-$r 1 r ). (5) 



*(roH^™. (6) 



To find the resistance we need to solve Eq. (||) . This is a Poisson-like equation and may be solved by using the lattice 
Green's function: 

V(r)=Rj2G(r-r')I(r>), (7) 

where the lattice Green's function is denned by 

A (r0 G(r-r') = -5 r ,r'. (8) 

Finally, the resistance between the origin and ro can be expressed by the lattice Green's function. Using Eq. (||) and 
(^) we obtain 

R(r ) = 2R[G(0)-G(r )], (9) 
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where we have made use of the fact that the lattice Green's function is even, i.e. G(r) = G(-r). Equation (|) is our 
central result for the resistance. 

To find the lattice Green's function defined by Eq. (||) we take periodic boundary conditions at the edges of the 
hypercube. Consider a hypercube with L lattice points along each side. Thus the total number of sites in the d 
dimensional hypercube is L d . Substituting the Fourier transform 



G w = ta £ G ( k ) e * kr ( 10 ) 

of the lattice Green's function into Eq. (@), we find 



L d 

kGBZ 



G(k) = i l y , (11) 
for the d dimensional hypercube where we have defined 

d 

e(k) = 2 (1 - cos ka,) . (12) 



Owing to the periodic boundary conditions, the wave vector k in Eq. ( |11[ ) is limited to the first Brillouin zone and is 
given by 

mi, 7719, rrid . 

k = -^b! + -^b 2 + • • • + -^b d , (13) 

where mi, m2, • • • rrid are integers such that —L/2 < rrii < L/2 for i = 1, 2, • • • , d, and by are the reciprocal lattice 
vectors defined by a;bj = 2nbi^ i, j = 1, 2, • • • , d. Here we assumed that L is an even integer, which will be irrelevant 
in the limit L — > 00. The mathematical desp^iptiaa-of the crystal lattice and the concept of the Brillouin zone can be 
found in many books on solid state physicstJiISEJO. 
Finally, the lattice Green's function takes the form 

kGBZ v 1 

If we take the limit L — > 00 then the discrete summation over k can be substituted by an integra 

1 f d d k 

L kGBZ J ^ BZ ( 27r ) 

where vq = a d is the volume of the unit cell of the d dimensional hypercube. Thus the lattice Green's function is 

JkeBZ (2n) £l K J 

Using Eqs. (^) and (|l^) in d dimensions the resistance between the origin and lattice point r is 



R(v ) = 2Rv (17) 
JkeBZ (27r) 

The above result can be simplified if we specify the lattice point as ro = /iai + £2^2 + ■ ■ ■ + Id&d'- 

* dxi r dx d 1 - e *(^i+-+^) 



R(h,k,---M = R -rrt- Trf ^7— -• (18) 



2tt 2?r 



Ei=i(l- cosa; i) 



From this final expression of the resistance one can see that the resistance does not depend on the angles between 
the unit vectors ai , a 2 , • • • , a^. Physically this means that the hypercube can be deformed without the change of the 
resistance between any two lattice points. The resistance in topologically equivalent lattices is the same. For further 
references we also give the lattice Green's function for a d dimensional hypercube: 

r d Xl r dx d e^^+-+'^> 
&{h,h,---,id)- ir ' / 77 ^ d n (1 9 ) 

J-TT Z7T J-TT Z7T 2 2^j =1 (l — COSXij 
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A. Conducting medium, continuum limit 



The resistance in an infinite conducting medium can be obtained by taking the limit as the lattice constant a 
tends to zero in Eq. ((l?]). Denoting the electrical conductivity of the medium by a, the resistance of a d dimensional 
hypcrcube with lattice constant a, according to Ohm's law, is given by 

R=-^- (20) 
aa a z 

Using the approximation e(k) w k 2 a 2 for | = a — ► 0, Eq. ( p7[ ) can easily be reduced to 

, s 2 f d d k l-e lkr ° , . 

R ^) = ~ ——d p ■ (2D 

The same result for the resistance of a conducting medium is given in Chaikin's bookEl. 

B. Linear chain, d = 1 

Consider a linear chain of identical resistors R. The resistance between the origin and site I can be obtained by 
taking d = 1 in the general result given in Eq. (|l8|): 

R(l) = R / =- ■ (22) 

t /_ 7r 27r 1 — cos x y ' 

The integral can be evaluated by the method of residues0 and gives the following very simple result: 

R(l) = Rl. (23) 

This can be interpreted as the resistance of I resistances R in series. The current flows only between the two sites 
separated by a finite distance. The two semi-infinite segments of the chain do not affect the resistance. 

C. Square lattice, d — 2 

Using Eq. jig) , the resistance in two dimensions between the origin and r = Ziaj + l 2 &2 is 

I"* dxi r d,xo 1 — (=.^121+^22:2 

R{hM) = R ^ P (24) 

J_ v 2tt _/_ 7r 2tt 2 — cos x\ — cos X2 

R r Axy_ r dx2_ 1 - cos (hxi + l 2 x 2 ) 

J_ n 2tt J_ 7r 2ir 2 — cos x\ — cos x 2 

The resistance between two adjacent lattice sites can easily be obtained from the above expression without evaluating 
the integrals in Eq. (|24|) or (pq). Note that interchanging l\ and l 2 in Eq. ( B5T ) does not change the resistance, i.e., 
R(h,h) — R(h,h)- This is consistent with the symmetry of the lattice. Then, from Eq. fl25|), the sum of R(0, 1) and 
i?(l,0) yields 



R(o,i) + R(i,o)= R r p- rp=R. 



(26) 



ThuSpfi{0,l) = R(1,0) = R/2. The resistance between two adjacent lattice sites is R/2, which is a well known 
resuhtxti 

In general the integrals in Eq. (|24|) have to be evaluated numerically. It is shown -in Appendix |a| how one integral 
can be performed in Eq. (|24|). We found the same result as VenezianQ and AtkinsorJj: 

J Q 7r smh s 
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where 



cosh s = 2 — cos y. (28) 



It turns out that this expression is more stable numerically than Eq. (25). As an application of the above formula 



one can calculate the resistance between second nearest neighbors exactly. After some algebra one finds 

R(l, 1) = R r — ^-^yf = 1 R . (29) 
J » * ^(2 -cosy) 2 -! * 

The energy dependent lattice Green's function of the tight-binding Hamiltonian for a square lattice is given by 

G{E-MM) = f ^ f ^ <*»d*i +'*»*) (30) 

J-TT 27T J-tt 27T — COS X\ — COS X2 



[see Eq. (5.31) in Economou's booka]. This is a generalization of our Green's function by introducing a new variable 
E instead of the value 2 in the denominator in Eq. ( |19|) for d = 2. Note that a factor 2 appearing in the denominator 
of our Green's function in Eq. (|l9|) is missing in Eq. (|30|). This is related to the fact that in the Schrodinger equation 
the Laplacian is multiplied by a factor of 1/2 while in our case the Laplace equation is solved. Comparing Eqs. ( JTs| ) 
(for d = 2) and (|30| ) one can see that the resistance is 

R(h,h) = R [G(2; 0, 0) - G(2; h,l 2 )} . (31) 

Based on the equation for the Green's function Morita0 derived the recurrence formulas for the Green's function 
G{E; li, I2) for an infinite square lattice (see Eqs. (3.8) and (4.2)-(4.4) in Morita's pape£3). Applying Morita's results 
(with E = 2) to the resistance given in Eq. (|3l]) we obtained the following recurrence formulas for the resistance: 

4m 2m — 1 

Rim + 1, m+ 1) = R(m, m) — R(m — 1, m — 1), 

v ' 2m +1 y ' 2m +1 y ' 

R(m + 1, m) = 2R(m, m) ~ R(m, m — 1), 
R(m + 1,0) = 4R(m, 0) - R(m - 1, 0) - 2R(m, 1), 

R(m + l,p) = AR(m,p) - R(m - l,p) - R(m,p + 1) - R(m,p - 1) if < p < m. (32) 

We have seen that R(l, 0) = R/2 and R(l, 1) = 2R/n. Since we know the exact values of R(l, 0) and R(l, 1) (obviously 
R(0, 0) = 0) one can calculate the resistance exactly for arbitrary nodes by using the above given recurrence formulas. 
This way we obtained the same results as Atkinson et. alfl using Mathematica. The advantages of our recurrence 
relations are that they provide a new, very simple and effective topi to calculate the resistance between arbitrary 
nodes on a square lattice. We note that van der Pol and Bremmera also gave the exact values of the resistance for 
nearby points in a square lattice using a different approach. 

It is interesting to see the asymptotic form of the resistance for large values of l\ or/and ig< In Appendix p] we 



derive the asymptotic form of the lattice Green's function for a square lattice [see Eq. (|B9|)]. Inserting Eq. (B9) into 
the general result of the resistance given in Eq. (||) the asymptotic form of the resistance is 

R(h, l 2 ) = - (in J 11 + q + 7 + ^ V (33) 



where 7 = 0.5772 ... is the Euler-Mascheroni constantlll. The same result was obtained by Veneziani except that we 
got an exact value of the numerical constant in Eq. ( p3[ ) whereas it was numerically approximated in Venezian's paper. 
The resistance is logarithmically divergent for large values of l\ and 1%. A similar behavior waSriouad for conducting 
medium by Chaikin et. al. in their bookt3. In the theory of the Kosterlitz-Thouless-BerezinskiiE2nl3 phase transition 
of the screening of topological defects (-scprtices) the same asymptotic form of the Green's function (las given in Eq. 
(|f39|)) has been used for a square latticed. Finally, we note that Doyle and Snell showed in their booka (pp. 122-123) 
that the resistance goes to infinity as the separation between nodes tends to infinity but the asymptotic form was not 
derived. 
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D. Simple cubic lattice, d = 3 



In three dimensions the resistance between the origin and a lattice point i"o = Ziai + ^28-2 + ^38-3 can be obtained 
from Eq. @: 



R{li,h, h) — R 



= R 



* dxi 



1 



a il\x\ +il 2 X2+U3X3 



d,X2 dxz 

2ir J_ n 2tt J_ 7r 2ir 3 — cos x\ — cos X2 — cos X3 
dx\ [* dx2 [* dx3 1 — cos (hxi + I2X2 + ^32:3) 



2tt 



2tt 



27r 3 — cos x\ — cos X2 — cos X3 



(34) 
(35) 



Similarly to the case of a square lattice, the exact value of the resistance between two adjacent lattice sites can be 
obtained from the above expression. Clearly, from Eq. (|35| ) (and for symmetry reasons, too) i?(l,0,0) = i?(0, 1,0) = 
i?(0,0,l) and 



i?(l,0,0) + i?(0,l,0) + i?(0,0,l) = i?r^i f ^ f p- = R. 

J-* 27r J-tt 27T J_ 7T 2n 



(36) 



Therefore, the resistance between adjacent sites is R/3. It is easy to show much in the same way that in a d dimensional 
hypercube the resistance between adjacent sites is R/d. 

Unlike in a square lattice, in a simple cubic lattice the resistance does not diverge as the separation of the en- 
tering and exiting sites increases, but tends to a finite value. One can write for the resistance R(h,l2,h) = 
2 [G(0, 0, 0) — G(h, h, h)] where G(l\,l2,h) is given in Eq. ( p9| ) for d = 3. It is well known from the theory of 
Fourier series (Riemann's lemma) that linip^oo dx(p(x)cospx — > for any integrable function <p(x). Hence, 
G(h,h,h) — * (which indeed corresponds to the boundary condition of the Green's function at infinity) and thus, 
R(h,hrh) — > 2G(0,0, 0) when any of,the h,l2,h — > 00. The value of G(0, 0, 0) was evaluated for the first time by 
Watsoncil and subsequently by Joyceca in a closed form in terms of elliptic integrals. The following exact result was 
found 



2G(0, 0, 0) = ( - ) (l8 + 12V2 - lOv^ - 7V6) [K(fc )] 2 = 0.505462 



where k = (2 — v / 3)(v / 3 — y/2) and 



K(fc) 



77/2 



d6 



(37) 



(38) 



is thexomplete elliptic integral of the first-kind. It is worth mentioning that a simpler result was obtained by Glasser 
et. al.c3 (see also Doyle's and Snell's bookQ) who calculated the integrals in terms of gamma functions: 



2G(0,0,0) = 



>/3 — 1 
96tt 3 



r 2 (i/24)r 2 (ii/24). 



(39) 



Thus, the resistance in units of R for a simple cubic lattice tends to the finite value 0.505462. . . when the separation 
between the entering and exiting sites tends to infinity. 

For finite separations the formula for the resistance given in Eq. (|3^) is not suitable for numerical purposes since 
the integrals converge slowly with increasing the number of mesh points. Similarly to the case of a square lattice, 
we can use the energy dependent lattice Green's function of the tight-binding Hamiltonian defined for a simple cubic 
lattice as 



" dxi 
2^ 



<V(E:/i./ 2 ./,>= / -7- I / -T- 



dx 3 cos (hxi + l 2 x 2 + I3X3) 



2tt E — cos X\ — cos X2 — cos X3 



(40) 



[see Eq. (5.49) in Economou's bookH]. This is a generalization of our Green's function by introducing a new variable 
E instead of the value 3 in the denominator in Eq. (|l9|) for d = 3. Note that the reason for the missing factor of 2 in 
Eq. @ is the same as that explained in the previous section. Comparison of Eqs. ( |l9| ) (for d = 3) and ( f40| ) yields 



R{hMM) = R [G(3; 0, 0, 0) - G(3; h, h M)\ 



(41) 



6 



The analytic behavior of the lattice Green's function G(E; l\, h, h) has-,been extensively studied over the past three 
decades. Numerical values of G(E; l\, 1%, Z3) were given by Koster et. al.Bj and Wolfram et. al.,E£lby using the following 
representation of the Green's function for a simple cubic lattice: 



G(E > 3;Z,m,n) 



dte-^h{t)I m (t)I n {t), 



(42) 



where 7/ is the modified Bessel function of order I. However, it turns out that another representation of the lattice 
Green's functions is more suitable for numerical calculations. For 3 < E the following forra-pf the Green's function 
along an axis was given in terms of complete elliptic integrals of the first kind by HoriguchiO 



G{E; 1, 0,0) 



dx K(fc)fc cos Ix, 



(43) 



where 



k = 



E — cos x 



(44) 



and K(fc) is defined in Eq. d38|). In Fig. |] we plotted the numerical values of the resistance using Eq. ( Q ) for E = 3 
and 1< I < 100. 



0.5 



R(1,0,0) 



0.45 



0.4 - 



0.35 - 



0.3 1 1 1 1 1 1 1 1 1 1 1 

20 40 60 80 100 

1 

FIG. 1. The resistance R(l,0, 0) in units of R along an axis for a simple cubic lattice. 

It is seen from the figme that the resistance tends rapidly to its asymptotic value given in Eq. (|37]) . We would men- 
tion that Glasser et. al.Eil gave other useful integral representations of the lattice Green's function for the hypercubic 
lattice for arbitrary dimension d. 

It was shown by Joycec3 that the function r £?(i?; 0, 0, 0) can be expressed in the form of a product of two complete 
elliptic integrals of the first kind. Horiguchicj and later on MoritaED obtained recurrence relations for the function 
G(E; I1J2, h) of the simple cubic lattice for an arbitrary site I1J2, h in terms only of GJ.E; 0, 0, 0), G(E; 2, 0, 0) and 
G(E; 3, 0, 0). The last two Green's functions were expressed by Horiguchi and Moritao in closed form in terms of 
complete elliptic integrals. From these results one can calculate the Green's function for an arbitrary lattice point. 
As an example, one of thqjecurrence formulas we obtained from the recurrence formula of the Green's function (at 
E = 3) given by Horiguch£j is 

R(m, 1, 0) = -- R(m — 1,0,0) + | R{m, 0, 0) - - R{m + 1, 0, 0). (45) 

Together with other recurrence formulas given in Horiguchi's and Morita's paper one can find R{1\, I2, h) for arbitrary 
values of h,l2,h- The numerical values of the resistance for small l\, ?2, ^3 are listed in Table I. 
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h,hM R{h,h,l 3 )/R 

Tfifi V3 

1.1.0 0.395079 

1.1.1 0.418305 
2,0,0 0.419683 

2.1.0 0.433599 

2.1.1 0.441531 
2,2,0 0.449352 
3,0,0 0.450372 
3,1,0 0.454416 
4,0,0 0.464884J— .r 

00 0.505462Elt 



TABLE I. Numerical values of R(h,l2,h) in units of R for a simple cubic lattice. 
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Our results are in agreement with those given by AtkinsonH. 

It is worth mentioning some useful references for the lattice Green's function of other three dimensional lattices, 
such as body centered and face centered. The exact values of.£j(l; 0,0,0) as well as other exact results can be found 
in Joyce's papeiO for a body centered cubic lattice. MoritaEZl gave the recurrence relations of the lattice Green's 
function for a body centered cubic lattice. InoueEJ derived the recurrence relations for a face centered cubic lattice 
and obtained the exact values of the lattice Green's function at some lattice points. Based on the results given in 
these papers one can calculate the resistance between arbitrary nodes for cubic lattices. 

Finally we would like to point out that in Eq. (|4l]) the resistance is related to the energy dependent Green's function 
at energy E = 3. JJpYever, it is known that the density of states at this energy corresponds to one type of the van 
Hove singularitiesEj'llj. Therefore, the resistance for a simple cubic lattice is related to the Green's function of the 
tight-binding Hamiltonian at the value of the energy at which the density of states has a van Hove type singularity. 
In the general case of a d dimensional hypercube the Green's function of the tight-binding Hamiltonian is infinite at 
E = d for d = 2 (logarithmically divergent) while it has a finite value for d = 3 as it was shown above. In the latter 
case the derivative of the Green's function with respect to E is singular. 



III. RECTANGULAR LATTICE 



In this section we shall calculate the resistance of a rectangular lattice in which the resistance of each edge is 
proportional to its length. Consider a rectangular lattice with unit vectors ai and a 2 and introduce parameter 
p = |ai| / |aa| . Let R be the resistance of the edge along the direction of a 2 . To find the resistance between the origin 
and site ro assume a current I enters at the origin and exits at site i"o. We denote the potential at lattice point r by 
V(r). Then, according to Ohm's and Kirchhoff's laws, we may write 



Al cct V(r) = -I(6 r , -6 r , ro )R, 



where the 'rectangular' Laplacian is defined by 



A T T ect V(r) 



V(r + ai) - V(r) V(r - ai) - V(r) 



V(r + a 2 ) - V(r) + V(r - a 2 ) - V(r). 



P P 

Then the equation for the lattice Green's function corresponding to the 'rectangular' Laplace equation is 

A™ ct G(r-r') = -5 r<r ,. 

The Green's function can be calculated in a way similar to the hypercubic case in Sec. 0. We have 

f d d k e 4kr 
Lebz {2v) d e(k) ' 



G(r) 



(46) 



(47) 



(48) 



(49) 



where vo — pa 2 is the area of the unit cell, and the Brillouin zone is a rectangle with sides 2tt/ |ai| and 2tt/ |&2 1 along 
the directions of a! and a 2 , respectively, and 



e(k) = 2 



(1 — cos kai) + 1 — cos ka 2 



Using Eq. (0), the resistance between the origin and lattice point r = li&i + ^ 2 a 2 for a given p is 



R(p;h,l 2 )=R 



R 



n dxi r dx 



_ n 2ir J_ 7T 2ir - (1 — cosa^) + 1 — cosa; 2 

71 dxi f 71 dx2 1 — cos (hxi + I2X2) 
2ir }_ v 2ir I (1 - cosxi) + 1 - cosx 2 



The integral over a; 2 can be evaluated similarly to a square lattice (see Appendix |a|) . We find 

dx l-e-^ s coshx 



R(p;h,l 2 )=R 



sinh 



(50) 

(51) 
(52) 

(53) 
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where 



coshs = H cosx. (54) 

P P 

Note that l\ and I2 are interchanged here compared to Eq. (|27|). 

Now it is interesting to see how the resistance between adjacent sites varies with p = |a.i | / |a.2 1 - The resistance 
between nearest neighbors is different along the x and the y axis. Thus, unlike in the case of a square lattice, 
R(p; 1, 0) 7^ R(p', 0, 1) for a rectangular lattice except for the trivial case p = 1 (square lattice). No sum rule exists for 
R(p; 1, 0) + R(p; 0, 1) as in the case of a square lattice. In Fig. || the resistances R(p; 1, 0) and R(p; 0, 1) are plotted as 
functions of p. 

1.7 r- 

1.5 - 
1.3 - 
1.1 - 



0.9 

0.7 - 
0.5 - 
0.3 - 

1 2 3 4 5 6 7 

P 

FIG. 2. The resistances R(p; 1, 0) and R(p; 0, 1) in units of R as functions of p for a rectangular lattice with p = |ai| / |&2| • 

From Fig. || one can see that resistance R{p] 1, 0) increases with increasing p. It can be shown that R(p; 1, 0) — » 00 
as p — > 00. This is physically clear since the lattice constant along the x axis increases resulting in an increasing 
resistance of each of the segments in this direction. On the other hand along the y axis a saturation of R(p; 0, 1) can 
be seen, which is not obvious at all. Expanding the integral in the expression of R{p] 0, 1) in powers of p, for large 
p we get R(p; 0, 1) « R(l - 2/np- 1 / 2 + 0{p~ 3 / 2 )). Thus, R(p; 0, 1) -> R when p — ► 00. Some additional results ax& 
presented for the energy dependent lattice Green's function for a rectangular lattice in the papers of Morita et. al.,L3 
and Katsura et. alEE 



IV. TRIANGULAR LATTICE 

In this section we calculate the resistance in a triangular lattice in which the resistance of each edge is identical, 
say R. First we consider a triangular lattice with unit vectors ai and a2, and with a lattice constant a = |ai| = |&2 1 - 
We choose ai and a.2 such that the angle between them is 120°. We introduce a third vector by a3 = — (ai +a2). The 
vectors drawn from each lattice point to its 6 nearest neighbors are ±ai,±a2,±a3. Assume that a current / enters 
at the origin and exits at site ro, and that the potential at site r is V(r). Again, from Ohm's and Kirchhoff's laws, 
we find 

A triang F(r) = _ j ^ _ ) R (55) 

where the 'triangular' Laplacian is defined by 

3 

&Y^V(r) = [V(r - a,) - 2V(r) + V(r + a*)] . (56) 

i=l 
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It is important to note that the triangular Laplacian A* riang defined above is 2/(3o 2 ) times that used for solving the 
Laplace equation with the finite-difference method on a triangular lattice. The factor 2/ (3a 2 ) ensures that the lattice 
Laplacian in the finite-difference method yields the correct form of the Laplacian in the continuum limit (a — > 0). 
However, in our case this limit does not exist since there are also connections (resistors) between adjacent lattice 
points along the direction of the vector a3 . 

The equation for the lattice Green's function corresponding to the 'triangular' Laplace equation is given by 

A*!' iang G(r - r') = -S r y. (57) 

The Green's function can be calculated in the same way as in the hypercubic case. We have 

f d 2 k p ikr 

G(r) = va 7^ rr» (58) 

JkeBZ (27r) 2 2^ i=1 (1 — cos ka,) 

where vo — \/3/2a 2 is the area of the unit celLand the vector k given in Eq. (|l3| ) runs over the Brillouin zone of the 
triangular lattice, which is a regular hexagor£3. 

The resistance between the origin and site ro can be obtained from Eq. (||) with the lattice Green's function for a 
triangular lattice and it yields 

„, . f d 2 k l-e lkr ° f d 2 k 1-coskro , , 

R(r )=Rv — ~2 3 — =Rv ~—2 ^ 3 7j j— ( 59 ) 

JkeBZ (27r) 2^ i=1 (1 — coska 4 ) JkeBZ (2ir) 2-ri=i U — coska;) 

Note that the factor 2 dropped out from the denominator. 

Using Eq. (|5^) we can easily find the resistance between two adjacent lattice sites. For symmetry reasons it is clear 
that R(&i) = R(sl2) = R(as). On the other hand 

Y,R&) = Rvof -^=R. (60) 



=i 



keBZ (27r)" 



In the last step we have made use of the fact that the volume of the Brillouin zone is 1/vo- Hence, the resistance 
between two adjacent lattice sites is i?(ai) = R/3. 

It is worth mentioning that without changing the resistance between two arbitrary lattice points one can transform 
the triangular lattice to a square lattice in which there are also resistors between the end points of one of the diagonals 
of each square. This topologically equivalent lattice is more suitable for evaluating the necessary integrals over the 
Brillouin zona-since the Brillouin zone becomes a square. The same transformation was used by Atkinson et. al. 
in their papertl. If we choose ai = (1,0) and ai = (0,1), then the resistance between the origin and lattice point 
ro = n&i + ma.2 is given by 

/*■ rl-r f rill 1 — P inx + im V 
/ 7T- o — ~ ? \ ' (61) 
_ n 2n 2n 3 — cos x — cos y — cos(x + y) 

In Appendix |c| we perform one integral in Eq. ( |6ll ) using the method of residuesS in a similar way to Appendix 
After some algebra it is easy to see that the result given in Eq. (C4) is in agreement with Atkinson's resuhi3: 

r' 2 dx l-e-'"- m '*cos(n + m)z 

R(n,m) = R / — , (62) 

j 7T smhs cos a; 

where cosh s = 2 sec x — cos x. 
For n — m in Eq. j62| ) we have 

C^ 2 dx 1 — cos 2nx 

R(n,n) = 2R ^^^^^^^^=. (63) 

Jo 77 ^/(3-cos2x) 2 -4cos 2 a; 

Evaluating this expression by means of the program Maple we obtained the same results for n = m as Atkinson!. 
Like in the case of a square lattice we believe that similar recurrence formulas exist for a triangular lattice but further 
work is necessary along this line. Similarly, it would be interesting to find the asymptotic form of the resistance as 
the separation between the two nodes tends to infinity. According to some preliminary work the resistance is again 
logarithmically divergent. 



11 



V. HONEYCOMB LATTICE 



In this section we calculate the resistances in an infinite honeycomb lattice of resistors. Atkinson and Steenwijki 
have studied this lattice structure exploiting the fact that the hexagonal lattice can be constructed from the triangular 
lattice by the application of the so-called A — Y transformation^! One advantage of our Green's function method is 
that we do not use such a A — Y transformation, specific only for the honeycomb lattice, therefore our method can 
also be used in a straightforward manner for other lattice structures. As we shall see later, in the honeycomb lattice 
each unit cell contains two lattice points, which is the main structural difference from the triangular lattice. Including 
more than one lattice point in the unit cell, the method outlined in this section can be viewed as a generalization of 
the Green's function method discussed in the previous sections. 

We assume that all the resistors have the same resistance R. The lattice structure and the unit cell are shown 
in Fig. ||. The angle between ai and a 2 is 120°, and |ai| = |&2 1 = VSa, where a is the length of the edges of the 
hexagons. There are two types of lattice points in each unit cell denoted by A and B in Fig. 0. 




FIG. 3. The honeycomb lattice with the unit cell, ai and &2 are the unit vectors of the lattice. Each unit cell contains two 
types of lattice points, A and B. 



From now on the position of the lattice points A and B will be given by the position of the unit cell in which they 
are located. Assume the origin is at one of the lattice points A, and then the position of a unit cell can be specified 
by the position vector r = na\ + ma.2, where n,m are arbitrary integers. We denote the potential and the current 
in one of the unit cells by V A (r) and Ve(r), and I A (r) and -/s(r), respectively, where subscripts A and B refer to 
the corresponding lattice points. Owing to Ohm's and Kirchhoff's laws the currents I A (r) and /s(r) in the unit cell 
specified by r satisfy the following equations: 



Ia(t) = 



V A {r)-V B {r) , Va(t)-Vb(t-bli) , V A (r) - V B (r - a x - a 2 ) 



R 



R 



R 



T V B {v)-V A {v) V b (t) - V A (r + ai ) V B (r) - V A (r + a x + a 2 ) 
Ib{v) = R + R + R ' 

Assuming periodic boundary conditions again the potential V^r) can be given by its Fourier transform 



(64) 



(65) 



keBZ 



where N is the-imrober of unit cells, and analogous expressions are valid for Vb(t), I A (r) and -Tb(i"). Here k is in the 
Brillouin zonetSEjO. Thus, we may rewrite Eq. (p3) as 



L(k) 



where 



V A (k) 
Ve(k) 



L(k) 



= -R 



IaM 



-3 h* 
h -3 



(66) 



(67) 
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and 



h= l+e lUai +e *k(ai+a 2 )_ (gg) 

Equation ( |66| ) is indeed the Fourier transform of the Poisson-like equation that determines the potentials for a given 
current distribution. However, in this case the Laplace operator is a 2x2 matrix. In r-represantion the analogous 
equation for a hypercubic lattice was given in Eq. (Q). 

The equation for the Fourier transform of the Green's function (which is a 2x2 matrix for a honeycomb lattice, too) 
can be defined analogously to a hypercubic lattice in Eq. (||) : 

L(k)G(k) = -1. (69) 

The solution of Eq. (|6^) for the Green's function is 

1 / 3 h* 

9-|fc| : 



G ( k ) = ^^U 3 )' 



where 9 — \h\ 2 = 2 (3 — coskai — coska2 — coskaa) and we have introduced a third vector, a3 = — (ai + &2)- 

There are four types of resistance. We denote the resistance between a lattice point A as the origin and site 
ro = nai +ma2 (which is an A-type site) by Raa(^o), while the resistance between the origin and site i"o + (2ai +a2)/3 
(which is a S-type site in the unit cell at ro) is denoted by Rab(?o)- For symmetry reasons, it follows that for the 
other two types of resistance: Rbb^o) = Raa(^o) and Rba(*o) = Rab {?(>)■ To measure Raa(^o) the current at sites 
A and B in the unit cell at r are 

I A (r) = I (6 r , - 5 T , ro ) and J B (r) = 0, (71) 
while for measuring Rab{to) we have 

I A (r) = IS rfi and I B {r) = -I5 r , ro . (72) 

Thus 

RMn) . K.(r = 0)-V,fa) . (73) 

RMta) _ v A ( r = o)-v a M (74) 



First, consider the resistance Raa{^o)- From Eqs. (p6|), (|69[)-(|7l|) we obtain 

y A (k)=/ J RG u (k)(l-e- ikr °). (75) 
Hence, using Eq. (|6|) we can obtain the r-dependence of the potential Vk(r) and then Eq. ( |73|) yields 



Raa(t ) = £ G "( k ) (! - coskr o) = *4 £ S , 1 C °t r ° ( 76 ) 

JV iV ' 3 — cos kai — cos ka 2 — cos ka 3 

keBZ keBZ 

Using Eq. (|l5|), the summation over k can be substituted by an integral and we obtain 

_ , . „„ f d 2 k 1 — coskr 

Raa(t ) = 3Rv / —2 r j- 2 r— , 77 

Agbz (27r) 3 - coskai - coska 2 - coska 3 

where vq = 3V3a 2 is the area of the unit cell. One can see that the same expression was found for a triangular lattice 
except for the factor 3«o here [see Eq. (|59|)]. Note that although the area of the unit cell vq is 6 times bigger than 
that for the triangular lattice, the area of the Brillouin zone is 6 times less and thus, the ratio of the resistances 
for honeycomb and triangular lattices is 3. Therefore, the resistance between two A-type nodes is three times the 
corresponding resistance in the triangular lattice. This again agrees with Atkinson's result. 

Similarly, using Eqs. (|6^), (p9|), ( |70| ) and ( [72] ) we can obtain VaC 4 ) and Vs(k). From Eq. ( B5[ ) the potential Va(v) 
can be determined, and analogously Vs(r). Finally, Eq. Mm leads to 
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R AB (ro) =RjjJ2 ( G n( k ) + G M - Gia(k) e^" - G 21 (k) e lkr °) 

kGBZ 




cPk 3 — coskro — cosk (ro + ai) — cosk (ro + ai + a 2 ) 
2n) 2 3 — cos kaj — cos ka 2 — cos ka 3 



= - [Raa (r ) + Raa (r + a x ) + Raa (r + ai + a 2 )] . (78) 

The same result was found by Atkinson and Steenwijk. 

The resistance between second nearest neighbor lattice sites may be found from Eq. ([77]). For symmetry reasons 
Raa(&i) — Raa{&2) = Raa(&3) and from Eq. ( f77j ) we have J2^=i Raa{sh) = 3i2. Thus, the resistance between 
second nearest neighbor lattice sites is R. Using Eq. (|78|), the resistance for adjacent lattice sites is Rab(0) = 
1/3 [Raa(0) + Raa(sli) + Raa(&i + a 2 )] = 2R/3, since obviously, Raa(0) = 0. We would just mention that the 
limiting value of -the resistance is again infinite as the separation between nodes tends to infinity (see p. 139 of Doyle's 
and Snell's book!). 

We note that one can transform the honeycomb lattice to a topologically equivalent brick-type lattice shown in 
Fig. without changing the resistance between two arbitrary lattice points. 



FIG. 4. The honeycomb lattice can be transformed into a topologically equivalent brick-type lattice. 
This fact has been utilized in the case of the triangular lattice in the previous section. 

As it can be seen from the above results, the resistance between arbitrary lattice points in a honeycomb lattice can 
be related to the corresponding triangular lattice. This kind of relation is the consequence of the so-called duality 
transformation in which the variables are transformed to Fourier transform variables (for more details see Chaikin's 
boolfiij pp. 578-584). The dual lattice of a triangular lattice is a honeycomb lattice. However, if the unit cell contains 
more than two non-equivalent lattice points then such a connection might not be used. On the other hand, the Green's 
function method outlined in the example of the honeycomb lattice can still be applied straightforwardly. 
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APPENDIX A: INTEGRATION IN THE EXPRESSION OF THE RESISTANCE FOR A SQUARE 

LATTICE 

Starting from Eq. ( |24| ) for the resistance i?(n, m) we have 

dy 



where 



R(n,m)=R[ ^ I(y), (AI) 



I(y) = I — 6 6 . (A2) 

_ 7r 2n 2 — cos x — cos y 
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Since R(n, m) — R(—n, m), we take n > 0. To proceed further, we perform the integral in I(y) using the method of 
residues. Introducing the complex variable z = e lx the integral I(y) reads 

/ ] 1 _ 7 n p irny 
7T n , n 5 , (A3) 
2tt 2z(2- cosy) - z 2 - 1 v ' 

where the path of integration is the unit circle. The denominator has roots at z + = e la+ and z_ = e la ~ , where a+ 
and a- satisfy the equation cos a — 2 — cosy and a_ = — a + . It is clear that for — n < y < tt we have 2 — cosy > 1, 
therefore the solution for a is purely imaginary. Thus we introduce s with a + = —a- = is, where s satisfies the 
equation 

cosh s = 2 — cos y. (A4) 

For — 7r < y < 7r it is true that s > 0, so the two poles of the integrand in I(y) are real numbers and satisfy the 
following inequalities: z + = e~ s <_L and 2_ = e s > 1. Thus the pole z + is within the unit circle, whil e Z- is outside. 
According to the residue theoremEj I(y) — —2i2iti ^2 residues within the unit circle. Using Eq. (A4) the residue of 
the integrand of I(y) at z + is 

1 l-e-™ s e imy 1 l- e - ns e mv , k , 

-, s = ■ (A5) 

27r2(2-cosy)-2z+ 2vr 2sinhs v ; 

We obtain 

1 1 _ „-ns imy i _ ~ns irny 

" ~ 2 ' 2 "2i 2 sinh 8 " sinh, ■ < A6 » 



where s satisfies Eq. (A4). Finally, the resistance R(n,m) for arbitrary integers n,m becomes 

r dy l-ejy» r dy 1 - e -H s cos my 

R{n,m) = R — = i? / — . (A7) 

J_ n Zir smns J ir smns 

The same result was found by VenezianB. 



APPENDIX B: THE ASYMPTOTIC FORM OF THE LATTICE GREEN'S FUNCTION FOR A SQUARE 

LATTICE 

In this Appendix we derive the asymptotic form of the lattice Green's function for a square lattice. The lattice 
Green's function at site r = is divergent since e(k) = for k = 0. Therefore we calculate the asymptotic form of 
G(0) — G(r). Starting from Eq. ( |l9| ) the lattice Green's function for site r = na.% + ma.2 in a square lattice becomes 

1 dv dr 1 - P inx P im v 

G(0) - G(n,m) = - — / — . (Bl) 

2 J_ 7r 2vr 7_ 7r 2vr 2 -cos x -cosy v ' 



The integral over x is the same as I{y) in Eq. (|A2j), so we can use the result obtained in Eq 

1 r dy l-e-WV" 1 " r dy 1 - e -H* cosmy 
G(0) - G(n,m) = - — — =/ — — . (B2) 



2 7-7T 2tt sinhs Jo 27r sinhs 

We follow the same method for deriving the asymptoticJjprm of the lattice Green's function for large m and n as 
Veneziantx A similar method was used in Chaikin's booktj (see pp. 295-296) in the case of a continuous medium in 



two dimensions. We break the integral in Eq. (|B2j) into three parts: 

G(0) - G{n, m) = h + h + h, (B3) 



where 



r7r dy 1-e-W* cosmy 

71 = ^ y ' 

"dy f 1 1 



o 2n \ sinhs y 
^ dy / e~\ n \ y cos my e~'"' s cosTOy 



o 2n \ y sinh; 



(B4) 
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and s satisfies Eq. (A4). The first integral I\ can be expressed by the integral exponential Ein(z) 



h = i Re 

Z7T 



dy 



I _ e (\n\-im)y 

y 



= -!-Re{Ein [it (\n\ — im)]} . 
2tt 



where Ein(z) is defined by 



Ein(z) = 



dt ■ 



t 



(B5) 



(B6) 



For large values of its argument, Ein(z) w lnz + 7, where 7 = 0.5772 ... is the Euler-Mascheroni constant. Thus, for 
large n and m I\ can be approximated by 



Ii ps — (In 1 7T (|n| — im)\ + 7) = — (in \J n 2 + m 2 + 7 + In 7^ 

Z7T Z7T V / 



Using Eq. (A4) the integral I2 can be evaluated exactly: 

1 




(2 — cos y) — 1 







)-i 









(B7) 



(B8) 



In the integral -Z3 the integrand is close to zero for small values of y and s since s ~ sinhs w y, while for larger 
values of y and s the exponentials are negligible, therefore 7 3 w 0. 

Finally, we find that the lattice Green's function for large arguments, i.e., |r| = a\Jn 2 + m 2 — > 00 becomes 



aw " - +^ 



(B9) 



The sameptesult is quoted on page 296 of Chaikin's bookEl. In the theory of the Kosterlitz-Thouless-Bprczinkii phase 
transitioiJi3ii3 Kosterlitz used the same asymptotic form of the Green's function for a square latticeO. 



APPENDIX C: INTEGRATION IN THE EXPRESSION OF THE RESISTANCE FOR A TRIANGULAR 

LATTICE 

Introducing the coordinate transformations x 1 = (x + J/)/2 and y' = (x — y)/2, and using the complex variable 
z = e ly 'I 2 we have R(n, m) — 1/2 R dx' /2tt I(x') from Eq. (|6l|), where 

— -, (CI) 

2tt z 2 cos ^- — z (3 — cos x') + cos ^ 

and the path of integration is the unit circle. The factor 1/2 in front of R(n,m) is the Jacobian corresponding to 
the transformations of variables. Since i?(nai — ma.2) = R(m&i — na.2), we take n — m > 0. The denominator has 
roots in z + = e la +l 2 and z_ = e ta -/ 2 and it is easy to see that they satisfy equations a + — — a_ and 2cos(a + /2) = 
(3 — cos x') I cos(ir'/2). For — 7r < x 1 < tt it is clear that a+ is purely imaginary. If we choose a+ such that Re a+ > 0, 
then z+ < 1 and z_ > 1. Thus, with a + = is, the residue of the integrand of I(x') at z + (which is inside the unit 
circle, i.e., s > 0) is 

1 1 _ -(n-m)f i(n+m)^- 11 — f> — (™ — m )l p i ( n + m )%' 

— ; = - — — , (C2) 

2ir 2z+ cos y — (3 — cos x 1 ) 2n 2 sinh § cos — 

where s satisfies the equation 

s 3 — cos a;' .„ , 

2cosh- = — . (C3) 

2 cos y 

Finally, the resistance for arbitrary integers n, m is 



R(n,m) = % ?L — ^ . (C4) 

2 7__ 27T sinh ^ cos 



16 



1 R. E. Aitchison, "Resistance between adjacent points of Liebman mesh," Am. J. Phys. 32, 566 (1964). 

2 B. van der Pol and H. Bremmer, Operational Calculus Based on the Two-Sided Laplace Integral (Cambridge University Press, 
England, 1955) 2nd ed., p. 372. 

3 P. G. Doyle and J. L. Snell, Random Walks and Electric Networks, (The Cams Mathematical Monograph, series 22, The 
Mathematical Association of America, USA, 1984) pp. 83-149. 

4 G. Venezian, "On the resistance between two points on a grid," Am. J. Phys. 62, 1000-1004 (1994). 

5 D. Atkinson and F. J. van Steenwijk, "Infinite resistive lattices," Am. J. Phys. 67, 486-492 (1999). 

6 After the bulk of this work was completed I received an e-mail from Sidney Redner who pointed out several references 
including Atkinson's and Steenwijk's workO. Reading these papers and demonstrating the usefulness of the Green's function 
method inspired me to include the recurrence formulas for an infinite square lattice and to show that our resistance formula 
for a triangular lattice agrees with that given by Atkinson and Steenwijk. 

7 S. Kirkpatrick, "Percolation and Conduction," Rev. Mod. Phys. 45, 574-588 (1973). 

8 E. N. Economou, Green's functions in Quantum Physics (Springer- Verlag, Berlin, Germany, 1983), 2nd ed., pp. 3-18, pp. 
71-96. 

9 S. Katsura, T. Morita, S. Inawashiro, T. Horiguchi, and Y. Abe, "Lattice Green's Function. Introduction," J. Math. Phys. 
12, 892-895 (1971). 

10 J. M. Kosterlitz and D. J. Thouless, "Ordering, metastability and phase transition in two-dimensional systems," J. Phys. C 
6, 1181-1203 (1973). 

11 J. M. Kosterlitz, "The critical properties of the two-dimensional xy model," J. Phys. C 7, 1046-1060 (1974). 

12 V. L. Berezinskii, "Destruction of long-range order in one- dimensional and two-dimensional systems having a continuous 
symmetry group I. Classical systems, "Sov. Phys. JETP 32, 493-500 (1971), [Zh. Eksp. Teor. Fiz. 59, 907-920 (1970)]. 

13 P. M. Chaikin and T. C. Lubensky, Principles of Condensed Matter Physics (Cambridge University Press, Cambridge, 
England, 1995), pp. 495-589, p. 46, pp. 100-101, p. 297, pp. 50-53. 

14 J. Lee and S. Teitel, "Phase transition in classical two-dimensional lattice Coulomb gases," Phys. Rev. B 46, 3247-3262 
(1992). 

15 J. M. Ziman, Principles of The Theory of Solids (Cambridge University Press, Cambridge, England, 1972), pp. 47-51, p. 12. 

16 N. W. Ashcroft and N. D. Mermin, Solid State Physics (Saunders College, Philadelphia, PA, 1976), p. 145 and p. 465, p. 89, 
p. 37. 

17 C. Kittel, Introduction to Solid State Physics (John Wiley and Sons, New York, 1986), 6th ed., p. 113, pp. 37-42. 

18 W. H. Press, B. P. Flannery, S. A. Teukolsky and W. T. Vetterling, Numerical Recipes (Cambridge University Press, 
Cambridge, England, 1990), pp. 615-622. 

19 G. B. Arfken and H. J. Weber, Mathematical Methods for Physicists (Academic Press, San Diego, CA, 1995), 4th ed., p. 421, 
p. 293, pp. 414-416. 

20 T. Morita, "Useful procedure for computing the lattice Green's function — square, tetragonal, and bcc lattices," J. Math. 
Phys. 12, 1744-1747 (1971). 

21 G. N. Watson, "Three triple integrals," Quart. J. Math. (Oxford) 10, 266-276 (1939). 

22 G. S. Joyce, "Lattice Green function for the simple cubic lattice," J. Phys. A 5, L65-68 (1972). 

23 M. L. Glasser and I. J. Zucker, "Extended Watson Integrals for Cubic Lattices," Proc. Natl. Acad. Sci. USA, 74 (1977), 
1800-1801. 

24 G. F. Koster and J. C. Slater, "Simplified Impurity Calculation," Phys. Rev. 96, 1208-1223 (1954). 

25 T. Wolfram and J. Callaway, "Spin-wave impurity states in ferromagnets," Phys. Rev. 130, 2207-2217 (1963). 

26 T. Horiguchi, "Lattice Green's function for the simple cubic lattice," J. Phys. Soc. Japan 30, 1261-1272 (1971). 

27 T. Morita, "Use of a recurrence formula in computing the lattice Green function," J. Phys. A 8, 478-489 (1975). 

28 T. Horiguchi and T. Morita, "Note on the lattice Green's function for the simple cubic lattice," J. Phys. C 8, L232-235 
(1975). 

29 G. S. Joyce, "Exact results for a body-centered cubic lattice Green's function with applications in lattice statistics. I," J. 
Math. Phys. 12, 1390-1414 (1971). 

30 M. Inoue, "Lattice Green's function for the face centered cubic lattice," J. Math. Phys. 15, 704-707 (1974). 

31 M. L. Glasser and E. Montaldi, "Staircase polygons and recurrent lattice walks," Phys. Rev. E 48, R2339-2342 (1993). 

32 T. Morita and T. Horiguchi, "Calculation of the lattice Green's function for the bcc, fee, and rectangular lattices," J. Math. 
Phys. 12, 986-992 (1971). 

33 S. Katsura and S. Inawashiro, "Lattice Green's function for the rectangular and square lattices at arbitrary points," J. Math. 
Phys. 12, 1622-1630 (1971). 

34 W. J. Duffin, Electricity and Magnetism (McGraw-Hill, London, 1990), 4th ed., p. 138. 

35 M. Abramowitz and I. Stegun, Handbook of Mathematical Functions (Dover Publication Inc., New York, NY, 1972), 9th ed., 



17 



p. 228. 



18 



